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Distribution of a two component electrolyte mixture between the model adsorbent and a bulk aqueous elec- 
trolyte solution was studied using the replica Ornstein-Zernike theory and the grand canonical Monte Carlo 
method. The electrolyte components were modelled to mimic the HCI/NaCI and HCI/CaCl2 mixtures, respec- 
tively. The matrix, invaded by the primitive model electrolyte mixture, was formed from monovalent negatively 
charged spherical obstacles. The solution was treated as a continuous dielectric with the properties of pure 
water. Comparison of the pair distribution functions (obtained by the two methods) between the various ionic 
species indicated a good agreement between the replica Ornstein-Zernike results and machine calculations. 
Among thermodynamic properties, the mean activity coefficient of the invaded electrolyte components was cal- 
culated. Simple model for the ion-exchange resin was proposed. The selectivity calculations yielded qualitative 
agreement with the following experimental observations: (i) selectivity increases with the increasing capacity 
of the adsorbent (matrix concentration), (ii) the adsorbent is more selective for the ion having higher charge 
density if its fraction in mixture is smaller. 
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1. Introduction 

A variety of phenomena of interest for basic science and technology occur when an electrolyte so- 
lution or a mixture of electrolytes is in equilibrium with a porous phase containing fixed charges. Com- 
prehension of the mechanism governing the equilibrium distribution of ions between the two phases is 
therefore of considerable interest. Materials containing charged micropores, such as membranes, glasses 
and certain clays, have a property of partially rejecting an electrolyte being filtered through them. In 
addition, such materials may be ion-selective; the ions having different sizes or charges adsorb to a dif- 
ferent amount. This makes a basis for the separation process used in water softening and purification, in 
treatment of radioactive wastes, and for analytical purposes as well. Ion-exchange is, therefore, not only 
a powerful tool in chemical analysis and in scientific research, but also of great importance for everyday 
life GUI. 

For a quantitative description of these phenomena, some knowledge of the spatial distribution of 
the counter-ion and co-ion species in the micropores is needed. To begin a mathematical description, 
a suitable model of the adsorbing material should first be assumed. The ion-exchange resins are quite 
often modelled as spherical holes or cylindrical channels. In such cases the calculations are, due to the 
known geometry, relatively simple (see, for example, 1 5]) and the Poisson-Boltzmann theory is most often 
accurate enough to make viable predictions of measurable properties. Alternatively, the Monte Carlo or 
Brownian Dynamics methods may be used for this purpose IbHlOll. 

In contrast to the models with well defined (spherical, cylindrical, etc.) geometry described above, 
we are here interested in modelling the porous systems where distribution of charges in a system is not 
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governed by a simple geometry. The porous material (matrix) is pictured as a set of charged obstacles. 
Distribution of the obstacles is obtained by a certain recipe (see, for example, |9j,|ll|]). Inhomogeneity 
that causes the field in which mobile ions are distributed is given on a molecular level. Such a model is 
a representative of a quenched-annealed systems which have been more extensively studied only in the 
last 25 years. The progress in this area of research is documented in several review papers I12l4l5ll . 

The pioneering contributions to the science of partly quenched systems are due to Madden, Glandt 
and coworkers (laJljl. Given and Stell HUS], Rosinberg 0], Bratko and Chakrab ortv ( isHiil. Chan- 
dler and coworkers |27|,|28|], Kaminsky and Monson 1 29, 30|], Holovko and coworkers l3lLl32fl . to mention 
only the most visible contributors. In references Il6lll711 and jioHi^l the replica Ornstein-Zernike theory 
was proposed to calculate the thermodynamic properties and spatial distribution of particles in partly 
quenched systems with short-range forces. 

Following these pioneering studies, Pizio and coworkers |[33j] succeeded in extending the replica the- 
ory to the systems with Coulomb interactions. In a series of papers liiUiill the quenched system was 
modelled as a set of charged obstacles containing positive and negative charges so that the matrix sub- 
system (the set of obstacles) was electroneutral. The other electroneutral subsystem, that is electrolyte 
solution, was assumed to occupy the void space in-between the obstacles, i.e., to invade the matrix. To 
calculate the thermodynamic properties and spatial distribution between the quenched and annealed 
particles the replica Ornstein-Zernike (ROZ) integral equations were applied. In addition, an expression 
was derived 1 38] for the mean activity coefficient valid within the replica hypernetted-chain approxi- 
mation. The canonical and grand canonical Monte Carlo method was used to verify the theoretical ap- 
proximations inherent to the replica theories. Very good agreement between the theory and simulations 
was obtained for both, thermodynamic parameters and for the pair distribution functions (see, for ex- 
ample, 1 37-39]). More recently, good agreement was confirmed by the independent Brownian Dynamics 
simulations based on the same model |9]. 

In the last years, the model for a quenched system was modified to allow studying the ion-exchange 
phenomena. The model proposed in references (ll|,|40j,|4l|] differs from the one presented above in one 
important aspect. We assume again that the distribution of charged obstacles corresponds to an equilib- 
rium distribution of ions in an electrolyte solution under conditions given by the temperature and the 
dielectric constant before quenching. However, in contrast to the previous studies, we assume that only 
anions are quenched and the cations are allowed to equilibrate (anneal). We also assume that the model 
system is in thermodynamic equilibrium with an external electrolyte solution, which penetrates into the 
matrix. A similar model, but for a somewhat simpler system, was previously studied by Pastore et al. [42]. 
The model was recently studied by the grand canonical Monte Carlo method |ll|,|4(j] as well as the Brow- 
nian Dynamics approach 11011 . The simulations indicate that the ROZ theory in the hypernetted-chain 
approximation, adapted to this model in Hill , yields very good agreement with computer simulations. 

In the present article, we extend the approach of the previous studies (ill Si t° the investigation 
of the ion-selectivity effect in the matrix representing a charged nanoporous system. The system of 
quenched obstacles formed from the primitive model ions is negatively charged. Within this system (ma- 
trix) the annealed ions are distributed; an excess of cations is present to maintain the electroneutrality 
condition. In equilibrium, the concentration of the invading electrolyte is determined by the properties 
of the bulk solution of the same chemical composition. Such a system can represent a crude model of the 
ion-exchange resin. The distribution of ionic species within the matrix is calculated by the ROZ integral 
equations and in several cases also by the grand canonical Monte Carlo (GCMC) method. For a certain 
composition of the external electrolyte, a different equilibrium composition of the internal electrolyte 
solution is established, favouring the ionic species, which interact stronger with the matrix particles. 
Structural, as well as the thermodynamic properties are calculated, such as the mean activity or Donnan 
exclusion coefficients, varying the composition of an external electrolyte. The ion-exchange isotherms, 
reflecting the ion selectivity in HCl/NaCl and HCl/CaCl2 mixtures, are presented as functions of the model 
parameters. 
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2. The model and theoretical methods 

To model an ion-exchange phenomenon, a quenched-annealed system composed of two subsystems 
was considered. The quenched negatively charged ionic obstacles represented the matrix subsystem, 
while the annealed subsystem was a mixture of two primitive model electrolytes with the addition of 
cations needed to compensate the negative charge of the matrix. It was assumed that the matrix remained 
insensitive to the presence of the annealed fluid, i.e., its structure did not change upon the introduction 
of the annealed phase. The solvent was treated as a structureless continuum (McMillan-Mayer level of 
description) characterized by the dielectric constant of pure water under the conditions of the study. 
Throughout the paper, indices and 1 designate the matrix and the annealed subsystems, respectively. 

The matrix phase was prepared in the following way: the +1 : -1 restricted primitive model elec- 
trolyte, with the diameters of ions o a + — a°_ — 4.25 A, was quenched at a certain temperature To, where 
the dielectric constant of the medium was £q. The structure of the matrix was considered to represent 
one of the equilibrium structures, corresponding to a given set of conditions (concentration, tempera- 
ture, relative permittivity), and governed by the Coulomb pair interaction potential 



( 1 j ' 1 j 

In equation \2.1\ , z9 (z°) is the charge of the ion i (j). equal to +1 or -1, Ab,o = e^/{47iE v Eok^To) is the 
Bjerrum length (eo denotes the unit charge, e v the permittivity of vacuum, and fce the Boltzmann con- 
stant), r designates the distance between the centres of ions i and j, and /3q = 1 / (fee To) . In the final stage 
of the matrix preparation, the cations were disregarded, so that the quenched phase consisted only of a 
spatially fixed arrangement of negatively charged (z°_ = — 1) particles. Conditions of the matrix prepara- 
tion corresponded to Ab,o = 7.14 A (water solution at 298 K), and ion concentrations cq = 0.1, 0.5, 1.0, and 
2.0 moldm" 3 . 

The annealed subsystem was modelled as a mixture of two + 1 : - 1 electrolytes, or a mixture of a + 1 : 
-1 and +2 : -1 electrolyte having a common anion. An additional number of univalent cations (common 
to cations of one of the electrolytes in the mixture), was present to neutralize the overall negative charge 
of the matrix. The diameters of the +1 ions were a\ = 5.04 A (model for hydrogen ions, H + ) or a\ = 3.87 A 
(model for sodium ions, Na + ), of the +2 ions a\ = 7.03 A (model for calcium ions, Ca 2+ ), and of the -1 
ions a\ = 3.63 A (model for chloride ions, CI"). Models of HC1 and NaCl, or HC1 and CaCl2 mixtures were, 
therefore, considered. Cations compensating the matrix charge were taken to be H + , mimicking in this 
way a H + -ion-exchange resins. The system was considered to thermally equilibrate under the conditions 
of Ti and ei, with the pair potential for a - 1 interaction being equal to 

and for a 1 - 1 interaction 

u f oo, r <(a)+o))l2, 
^ U H=\ z}z)X^lr, r>&rJ + 4/2. (2 " 3) 

In equations J2.2t and (2.3) , Ab,i = e\l ' {\ne v e\k^,T\] is the Bjerrum length for the annealed solution, and 
Pi - l/(A;Bri). Although the conditions for the preparation of the quenched phase can differ from the 
conditions of observation, we assumed them, to be equal in this study, i.e. Ag ( i — Aj^o — 7.14 A. 



r<{a° + <j°)/2, 



2.1. The replica Ornstein-Zernike integral equation theory 

A case of adsorption of a single electrolyte within the charged matrix was considered by our group 
several years ago 1 11]. Here, we extend the theoretical description to the study of mixtures of electrolytes 
with a common ion. 

In order to define the spatial distribution of matrix ions we first need to obtain the structure for a 
+ 1 : - 1 electrolyte. This is obtained by solving a set of Ornstein-Zernike (OZ) equations in the form 
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where h and c stand for the total and the direct correlation functions, respectively, and ® denotes the 
convolution in r-space. Due to spherical symmetry, h°P_ = h°P + and c°_°_ = c°_ ° + . Equations 12.4> need to be 
renormalized before being solved numerically (see, for example, The hypernetted-chain (HNC) 

closure relation between h and c (c,j = h{j -ln[/j i; - + l] - PoU??) was used to obtain the solution. Note 

that only is needed in the subsequent replica OZ equations 12.5\ . This function contains all the infor- 
mation of the matrix subsystem. 

The replica Ornstein-Zernike equations for the quenched-annealed system (cf. |22|,[4l|]) read 

1,10 „io . 10 _ „o u oo , v r^n ^ „i 1,10 „i2_„i 1,10 1 ,r, c v 

h i _=c l _ + c i _®p_h__+ \c ik ®p k h k _-c ik ®p k h k _\, (2.5) 

fc=A,B,- 

h\) = eg + c» . p°_ ft!* + J£ J C U 8 p l h ll _ C 12 . p l fc 21 ] , (2 . 6) 

tf = c% + c?_®p°_h°) + £ [^•pihg + cg.pifcg-Z^.p^], (2.7) 

fc=A,B ,- 

where indices i and j stand for cations or anions of the annealed electrolytes in the mixture (i.e., H + , Na + 
and Cl~, or H + , Ca 2+ , and CD, while A and B denote only the cations of the annealed mixture (H + and 
Na + , or H + and Ca 2+ , respectively). We stress that the pjj + is the total density of H + ions in the system 
which belong to the HQ and to the cations compensating the negative charge of the matrix. 

As in the case of matrix preparation (see expression H2.4» . equations 12.5M2.7> need to be renor- 
malized prior to numerical solution (i.e. c and h need to be split into a short and long range part). The 
procedure is documented in (34j[3S|]. The only difference from the published equations is that here k 2 
contains the sum over all the annealed species. Additionally, HNC closure condition was applied to corre- 
late the h and c functions: 

c mn = h nin_ ln \ h mn + jl _ fl ^m* (2 g) 

c}] = h}]-\n\h} 2 j+ l\, (2.9) 

where m and n take the values and 1, and i (j) in equation )2.8( corresponds to the annealed cation 
or anions (H + , Na + , Ca 2+ , and Cl~) and to the matrix ions, i (j) indexes in equation <2.9) belong only 
to annealed ions. Note that ujf = since there is no interaction between the ions of different replicas. 
Equations <2.5Hl2"77l with the HNC closure were solved by a direct iteration on a grid of 2 15 points with 
Ar = 0.005 A. 

The individual activity coefficient within the replica HNC formalism reads (i^] 

^T) - - P °-[ cJoVdr- E P)f[^i-^;]dr 

J ;'=a,b- J 

+ y { ^ k - a dr 4 j } B _ p ) ! ( *s H - c ") - h u ( h u - c ") 1 dr - (2 - io) 

where q S ) denotes the short range part of the direct correlation function coming from the renormal- 
ization procedure, and dr = 4nr 2 dr. Again, A and B denote cations in the mixture (H + and Na + , or H + 
and Ca 2+ ). 

2.2. Monte Carlo simulations 

The matrix was prepared from a size symmetric (cr° = cr°_ = 4.25 A) +1 : -1 electrolyte using a sepa- 
rate canonical Monte Carlo simulation. After equilibration, the matrix anions were frozen in their posi- 
tions, while the cations were replaced by the model hydrogen cations (<7h+ = 4.25 A) and made it possi- 
ble to move freely, equilibrating together with the annealed electrolyte mixture (HCl/NaCl or HCl/CaCl2). 
The annealed electrolyte ions were then distributed within the matrix and the system was studied us- 
ing the Monte Carlo method in the grand canonical ensemble. The methodology of this approach is well 
established and extensively described in several previous papers and is not, therefore, repeated here 
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The details of simulations are as follows: The number of matrix particles in the simulation box varied 
from 300 to 2500, depending on the concentration. The average number of fluid cation species distributed 
within the matrix varied from 500 to 3500. The periodic boundary conditions within the minimum image 
convention were applied. The ions within the matrix were first equilibrated over at least 10 7 Monte Carlo 
steps. After the equilibration, a production run of 2 • 10 8 attempted configurations was carried out to ob- 
tain the average concentration of the adsorbed electrolyte species. Note that thermodynamic properties 
were averaged over different annealed fluids, as well as over one to two different matrix configurations. 

To calculate the equilibrium concentration in the bulk electrolyte, the activity coefficients of the bulk 
electrolyte mixture should be known in advance. These values were obtained using the hypernetted- 
chain (HNC) theory which has proved to be very successful in describing the properties of ionic fluids |45]. 

3. Results and discussion 

The two theoretical methods described above were established to be complementary in the sense 
that they yield consistent results for structural and thermodynamic properties of electrolyte solutions 
adsorbed in nanoporous materials (n[[3^[38l[4(l[4l]]. While the ROZ/HNC theory is computationally less 
demanding, and is, therefore, faster to use for systematic investigations of different effects on the system 
properties, the GCMC method is more convenient in studying the equilibrium distribution of ions between 
the bulk and matrix phase. In this paper we used both methods; the results are organized as follows. First, 
we tested the performance of the ROZ/HNC theory for the case studied here by comparing the structural 
(pair distribution functions) and thermodynamic properties (mean activity coefficients, y+ + v ~ = y + + j-~ ) 
obtained using the two methods. Next, the effect of the matrix concentration (e.g. adsorbent capacity) on 
the mean activity coefficient, which determines the chemical equilibrium in the system, was examined 

Table 1. The mean activity coefficients of electrolytes adsorbed in the matrix: y+, obtained from the 
ROZ/HNC theory and GCMC simulation for different matrix concentrations. Concentrations of HC1 and 
NaCl (chci and c Na Q, respectively) and of the matrix (co) are given in moldm -3 . 7 0Ut = 0.5 moldm -3 , 
A B = Ab,i = 7 -14 A. The numerical error of both methods was estimated to be in the last digit given. 
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0.903 


0.930 




0.0325 


0.202 


1.057 


1.086 


0.906 


0.932 




0.0206 


0.240 


1.064 


1.092 
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0.874 


0.883 
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0.798 




0.0562 


0.277 


0.872 


0.880 


0.789 


0.797 




0.0325 


0.321 


0.874 


0.881 


0.789 


0.796 




0.396 


0.0494 


0.805 


0.806 


0.755 


0.754 




0.347 


0.0990 


0.801 


0.802 


0.750 


0.750 




0.298 


0.148 
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0.396 
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0.774 
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using the ROZ/HNC theory. The Donnan exclusion coefficients and the ion-exchange isotherms, reflecting 
the selectivity of an adsorbent were evaluated by the Monte Carlo simulation in the grand canonical 
ensemble. 



3.1. Comparison of the ROZ/HNC results with computer simulations 

To test the accuracy of the ROZ/HNC theory for the model of interest, we compared the pair distri- 
bution functions and the mean activity coefficients with the results of computer simulations. Figures [l] 
and[2]show the pair distribution functions (++ and H — ) of the annealed ions, while figure [3] shows the 
annealed ion-matrix distribution functions and -J€). In all cases the concentration of the matrix 
was Co = 1.0 mol dm -3 . Figure[T]and the left panel of figure[3]apply to the mixture of HC1 and NaCl, while 
figure [2] and the right panel of figure [3] describe the distributions in the mixture of HC1 and CaCl2. Sym- 
bols denote the GCMC and the lines denote the ROZ results. Excellent agreement of the two methods is 
obtained for distribution functions of hydrogen, while the agreement is slightly worse, but still good, for 
the functions belonging to sodium and calcium ions. Note that the hydrogen ions are the ones compen- 
sating for the adsorbent charge and are, therefore, present in an excess to other ions. A similarly good 
agreement is obtained for the annealed ion-adsorbent pair distribution functions shown in figure[3] 




5 10 15 20 5 10 15 20 

r/A r/A 

Figure 1. Fluid-fluid pair distribution functions (mixtures of HC1 and NaCl). Left: H + -H + (•) and H + - 
Cl" (O); Right: Na + -Na + (•) and Na + -Cl" (O). c = 1.0 mol dm" 3 , cgL = 0.0325 mol dm" 3 , cjj a = 
0.2018 mol dm" 3 . Symbols represent GCMC data while lines show the ROZ/HNC theory. A B = A B>1 = 
7.14 A. 




5 10 15 20 5 10 15 20 

r/A r/A 

Figure 2. Fluid-fluid pair distribution functions (mixtures of HC1 and CaCFi). Left: H + -H + (•) and H + - 
Cl" (O); Right: Ca 2+ -Ca 2+ (•) and Ca 2+ -Cl" (O). c = 1.0 mol dm" 3 , c™ a = 0.0436 mol dm" 3 , c™ a = 

0.0522 mol dm" 3 . Symbols represent GCMC data while lines show ROZ/HNC theory. A Bj0 = A B>1 = 7.14 A. 
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5 10 15 20 5 10 15 20 

r/A r/A 

Figure 3. Fluid-matrix pair distribution functions (mixtures of HC1 and NaCl (left) or HC1 and CaCl2 
(right)). Left: H+-./7 (•), Na + -.<7 (O), and CT-M (■), c£L = 0.0325 moldm" 3 , c™ aC{ = 0.2018 moldm" 3 ; 
Right: Y^-Ji (•), Ca 2+ -^ (O), and QT-Ji (■), c™ Q = 0.0436 moldm" 3 , = 0.0522 moldm" 3 . 

cq = 1.0 mol dm" 3 . Symbols represent GCMC data while lines show ROZ/HNC theory. Ab,o = <^b,1 = 7 -14 A. 



Table 2. The same as in table[TJ but for HCl/CaCl2 mixtures. 
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0.808 
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0.0952 
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0.00804 


0.0880 
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1.0 
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1.047 


1.074 
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0.801 




0.0253 


0.0675 


1.043 
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0.765 
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0.0436 


0.0522 
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0.0751 


0.0308 
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0.00221 
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0.621 


0.624 




0.00585 
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0.0403 
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0.860 


0.867 


0.614 
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0.0670 
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0.861 


0.869 


0.617 


0.627 




0.108 
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0.865 
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0.620 


0.631 
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0.628 


0.638 




0.00296 
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0.542 
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0.0126 
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0.763 


0.544 
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0.0291 
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0.765 


0.765 


0.545 


0.546 




0.0524 


0.135 


0.767 


0.767 


0.547 


0.548 
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0.0827 


0.123 
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0.770 


0.550 


0.551 




0.121 


0.120 


0.774 


0.775 


0.554 


0.555 




0.171 


0.093 


0.780 


0.780 


0.558 


0.559 




0.236 


0.071 


0.787 


0.787 


0.564 


0.565 




0.323 


0.0411 


0.796 


0.797 


0.572 


0.572 



The results for the mean activity coefficients, j±, in aqueous mixtures of HQ with NaCl, and HQ 
with CaCl2 in the matrix are collected in tables [T] and [2] respectively. Note that in all cases, the partic- 
ular electrolyte concentrations given, are in equilibrium with the bulk solution of ionic strength / out = 

0.5£ ; z?c out = 0.5 moldm" 3 . The agreement between the ROZ/HNC and GCMC results is excellent for 

1 — ^ 

small matrix concentration, but it deteriorates for larger en values. At cq = 1.0 moldm , the mean activ- 
ity coefficients calculated within the ROZ/HNC theory are a bit too large, as already noticed in previous 
papers (ll],[36|,[38|,|4(J. Nevertheless, the differences discussed here are within the estimated numerical 
errors of the two methods. 
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3.2. The effect of the adsorbent capacity on the mean activity coefficient of the invad- 
ing electrolyte. 

The ROZ/HNC method was used to investigate the effect of the adsorbent capacity, given here as the 
matrix concentration, on the mean activity coefficient of the invading electrolyte. Note that while the 
ionic strength of the annealed electrolyte mixture within the matrix (denoted by superscript "in") was 
kept constant, I m = 0.5 moldm -3 , for all cases examined here, the composition given by 



Y m — 

a hci _ 



HCl 



L HC1 L Z 



= 1 



(3.1) 



varied. Z in equation 13. 1> denotes either NaCl or CaCl2. Concentrations of the matrix obstacles considered 
were: c = 0.1, 0.5, 1.0, and 2.0 moldm" 3 . The ROZ/HNC results are shown in figures H(HCl/NaCl), and[5] 
(HCl/CaQ.2) by continuous lines, the dashed lines (OZ/HNC) show the mean activity coefficients in the bulk 
solution with I = 0.5 mol dm -3 . 
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Figure 4. The mean activity coefficient as a function of the (mixtures of HCl and NaCl), obtained 
using the ROZ/HNC theory. Left: HCl; Right: NaCl. From top to bottom (continuous lines): c = 2.0, 1.0, 0.5, 
and 0.1 moldm -3 . Dashed line (obtained using OZ/HNC) represents the bulk value. 7 m = 0.5 moldm -3 , 
A B ,0 = ^B,1 =7.14 A. 
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Figure 5. The same as in figure|4] but for HCl/CaCl2 mixtures. 



As the matrix concentration (capacity of the adsorbent) increases, the mean activity coefficients of 
all the species in the solutions increase from values smaller than unity to values higher than unity, as 
previously observed for a single electrolyte adsorbed in charged matrix (lll|40j]. The matrix charge re- 
pels the co-ions of the invading electrolyte mixture, which causes the exclusion of the electrolyte on the 
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Figure 6. Donnan exclusion coefficient as a function of X™ a (mixtures of HC1 and NaCl), obtained with the 
GCMC simulations. Left: HC1; Right: NaCl. c = 1.0 (•), 0.5 (O), and 0.1 (■) moldm" 3 . / out = 0.5 moldm" 3 , 
Ab,0 = %1 =7.14 A. 



whole from the matrix. Concentration of a particular electrolyte is smaller in the matrix than in the cor- 
responding equilibrium bulk solution. This can be clearly seen from the values of the Donnan exclusion 
coefficient, T,-, for a particular salt i (HQ, NaCl, or CaCl2). Tj is defined as 



(3.2) 



The results for the Donnan exclusion coefficients, calculated by the GCMC method (ionic strength of 
the bulk electrolyte was kept constant, 7 0ut = 0.5 moldm -3 ), are shown in figures [6] (HCl/NaCl) and [7] 
(HCl/CaCl2). As indicated by the behaviour of the mean activity coefficients (figures|4]and[5}, the annealed 
electrolytes are excluded from the adsorbent (positive T values) in all cases. The exclusion is stronger 
for higher matrix concentration (i.e. higher adsorbent capacity, or higher charge), and depends on the 
composition of the invading electrolyte. The exclusion coefficient for all electrolytes in the mixture de- 
creases with the increasing fraction of the HC1 in the solution. While the r H ci approaches the value one 
as -X^jL — ► (i.e. for no exclusion), the exclusion coefficient of the other electrolyte in the mixture ap- 
proaches a somewhat lower value. Such a result has previously been obtained for electroneutral adsor- 
bents |41]. The difference in the behaviour of the two electrolyte components is caused by the modelling 
of the system (see section 2), mimicking the H + -ion-exchange resin. The matrix charge is namely "neu- 
tralized" with hydrogen ions, which are accordingly always present in the annealed electrolyte mixture. 





Figure 7. The same as in figure[6] but for mixtures of HC1 and CaCl2. 
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3.3. The ion-exchange isotherms 

Due to a different degree of adsorption (exclusion), the concentration ratio of the competing ions in 
the adsorbent (cations) is generally different from that in the bulk solution. In other words, the adsorbent 
"selects" one species in preference to the other. The study of selectivity was one of the main goals of this 
work. Here we present the results in the form of ion-exchange isotherms, showing the mole fraction of 
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Figure 8. Ion-exchange isotherms (mixtures of HC1 and NaCl), obtained with the GCMC simulations. Left: 
HC1; Right: NaCl. c = 1.0 (•), 0.5 (O), and 0.1 (■) moldm" 3 . A B ,o = %i = 7.14 A. Lines are guides for the 
eye. 



one electrolyte component within the adsorbent, X m , as a function of the composition of the equilibrium 
external solution, X out . The ion-exchange isotherms are presented in figure [8] (HCl/NaCl mixture) and 
figure[9](HCl/CaCl2 mixture). Symbols show the grand canonical MC results at different matrix concentra- 
tions, and the connecting lines merely serve to guide the eye. The diagonal lines apply to a hypothetical 
case with no preference to the adsorbent. 





Figure 9. The same as in flgure|8] but for mixtures of HC1 and CaCl2. 



As noticed before |41, 46], the selectivity increases with the increased capacity of the adsorbent (in- 
creased matrix concentration). In all cases, the mole fraction of the HQ solution in the adsorbent is 
smaller than in the external solution (left panels of figures [8] and [9), the opposite is true for the NaCl 
or CaCl2 (right panels of figures [8] and [9), respectively. The selectivity strongly depends on the composi- 
tion of the solution as well. For the component that is preferentially adsorbed in the matrix (NaCl and 
CaCi2 in our case), the selectivity increases with an increasing fraction of the component in the solution, 
while the opposite is true for the other component. Similar trends were observed experimentally fUd]. 
Due to a stronger electrostatic interaction between the matrix charges and the ion with higher charge 
density (Na + or Ca 2+ , respectively), the latter ions are excluded from the matrix to a smaller extent. 
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4. Conclusions 

The model of ion-exchange resin was studied using the ROZ/HNC theory and Monte Carlo simulation 
in grand canonical ensemble. Theoretical results for the pair distribution functions were found to be in 
good agreement with computer simulation results; it is confirmed that ROZ theory represents a viable 
alternative to computer simulations. For all the electrolyte mixtures and matrix concentrations studied 
here, the electrolyte components are excluded from the adsorbent phase. This is most often explained as 
a consequence of electrostatic repulsion between matrix charges and co-ions from the electrolyte com- 
ponent. However, the ionic species are adsorbed to a different amount. Consistently with experimental 
observations, the ions with higher charge density (sodium and calcium ions in our case) adsorb to a 
greater extent. 
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Mofle/iiOBaHHfl ioHHO-o6MiHHoT B3aeMOfliT b Haiionopi/iCTi/ix 
MaTepia/iax 

M. JlyKiuiH, B. B/iaxi/i, B. Xpi6ap-/li 

YHiBepcuTeT /lKj6/iflHW, <|)aKy/ibTeT xiMii' Ta xiMiHHi/ix TexHO/ioriPi, AuiKepneBa 5, SI-1000, f\K>6mHa, OiOBeHifl 
BuKopucTOByKDHM Teopiro perviiKOBaHoro piBHflHHH OpHLUTetiHa-l4epHiKe i MeTOfl MoHTe Kapyio y Be/inKO- 

My KaHOHiHHOMy aHCaM6fli, flOOlifl>KyETbCfl p03nOfli/l flBOKOMnOHeHTHOi e/ieKTpO/liTMHHOI CyMil±li Mi>K MOfle/lb- 

hi/im aflCop6eHTOM i BOflHMM po34HHOM e/ieKTpo/iiTy. KoMnoHeHTH eyieKTpoyiiTy MOfle/iKJKJTb cyMiwi HCI/NaCI i 
HCI/CaCl2. MaTpni4fl, 3anoBHeHa npuMiTHBHOK) Mop,ennK) e/ieKTpo/iiTHMHoi cyrviiiui, 6y/ia ccf>opMOBai-ia 3 moho- 
Ba/ieHTHHx HeraTHBHO 3apfl,q>KeHHx ccf>epn4Hi<ix HacTHHOK-nepeinKOfl. Po34hh po3r/iflflaBCfl ak HenepepBHwfi 
flie^eKTpnK 3 B/iacn/iBOCTflMn 4wctoi' BOfln. nopiBHAHHfl napi-tux q>yHKi4iM po3nofli/iy (oTpi/iMaHe o6oMa Me- 
TOflaMn) Mi>K pi3HHMH iOHHMMH copiaMW BKa3yK)Tb Ha flo6pe y3roflxeHHfl Mix pe3yyibTaTaMH Teopii' pen/iiKi/i 
OpHi±iTeMHa-l4epHiKe i MaiunHHwx po3paxyHKiB. CepeflTepMOflUHaMiHHux B/iacTHBOCTeii, po3paxoBaHo Koecjsi- 
LiieHT cepeflHboT aKTHBHOCTi KOMnoHeHTiB e;ieKTpo;iiTy b MaTpnuj. 3anponoHOBaHO npody MOfle/ib pj\n iOHHO- 
o6MiHHoT B3aeMOfliT. Po3paxyHOK ce/ieKTHBHOCTi flaB nKiCHey3rofl>KeHH5i 3 TaKHMn eKcnepuMeHTa^bHUMw flaHw- 
mu: (i) ce/ieKTi/iBHicTb 3pocTae 3 poctom CMHOCTi aflcop6eHTa (KOHu,eHTpau,ii Maipnuj), (ii) aflcop6eHT npoHB/iae 
6mbiuy ce/ieKTHBHiCTb p,o ioHy 3 buluokj 3apaflOBOK) rycTHHOK), AKU40 fioro nacTKa b cyMiiui e MemuoK). 

K/iKJMOBi c/ioBa: ap,cop6u,tn, 3MimaHi ejieKTpojiim, ceneKmBHiCTb, MoHTe Kapno, reopw penniKM 
OpHi±iTetit-ia-l4epHiKe 
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